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We derive the constitutive equations of causal relativistic dissipative hydrodynamics [d- 
hydrodynamics) from perfect nonextensive hydrodynamics (g-hydrodynamics) using the nonexten- 
sive/dissipative correspondence (NexDC) proposed by us recently. The g-hydrodynamics can be 
thus regarded as a possible model for the d-hydrodynamics facilitating its application to high en- 
ergy multiparticle production processes. As an example we have shown that applying the NexDC 
to the perfect 1 + 1 g-hydrodynamics, one obtains a proper time evolution of the bulk pressure and 
the Reynolds number. 

PACS numbers: 24.10.Nz, 25.75.-q 



The ideal hydrodynamic model successfully reproduces 
most of the RHIC experimental data [1,] suggesting that 
matter created there resembles a strongly interacting 
quark-gluon fluid Q rather than a free parton gas. How- 
ever, there are also hints that, after all, a hadronic 
fluid cannot be totally ideal Q and should be described 
by some kind of dissipative hydrodynamic model (or 
d- hydrodynamics) 0, 0, 0]- A number of such mod- 
els have been proposed recently 0, i i [13, Eii. Al- 
though very promising, they have many difficulties, both 
in what concerns their proper formulation and appli- 
cations [1, [13, Il2|- Any new attempt to understand 
them more deeply is therefore welcome. In this letter we 
demonstrate that one can obtain constitutive equations 
of causal d-hydrodynamics starting from the relativistic 
perfect nonextensive hydrodynamics (q-hydrodynamics) 
supplemented by the nonextensive/dissipative correspon- 
dence (NexDC) proposed in [l^ and taking place be- 
tween the perfect q-hydrodynamics (based on the nonex- 
tensive Tsallis statistics |14] and described by the nonex- 
tensivity parameter q) and the usual d-hydrodynamics 
(based on the extensive Boltzmann-Gibbs statistics (BG) 
to which Tsallis statistics converges for 5 1). We 
argue that the d-hydrodynamics emerging from such q- 
hydrodynamics is causal and preserves the simplicity of 
the latter which facilitates numerical applications [l3| 
(instead of solving very complicated second order equa- 
tions of d-hydrodynamics one can solve the much simpler 
equation of motion of the q- hydrodynamics). 

Our q-hydrodynamics [isj is based on the relativis- 



tic nonextensive kinetic theory proposed in [15| (see also 
[1^). It is understood that such a theory accounts, by 
means of the parameter q, for all kinds of possible strong 
intrinsic fluctuations and long-range correlations existing 
in a hadronizing system. It replaces the usual notion of 
local thermal equilibrium by a kind of stationary state 
which also includes some interactions. It is described by 
the Tsallis nonextensive statistics [3l and characterized 



by the nonextensive parameter q > 1 [17|, |18|, ll9|. In 
this approach the perfect q-hydrodynamic equations (i.e., 
equations for the perfect nonextensive g-fluid without any 
additional currents) are given by the g- version of the en- 
ergy momentum tensor: 



= 0. 



(1) 



Here eq{Tq), Pq{Tq), Tq = Tq{x) and u'^{x) are, respec- 
tively, the nonextensive energy density, pressure and tem- 
perature field [isj and an accompanying hydrodynamic 
q-flow four vector field (which allows decomposition ^ 
to be performed; A^'' = g^''' - u^^u^q). If T{x) and M^(a;) 
are the corresponding temperature and velocity fields in 
the case of BG statistics {q — 1) with e and P being the 
usual energy density and pressure then, as proposed in 
[l^ . one can map a q-Aow into some dissipative d-flow 
by requiring that the following relations hold: 



3n, 



(2) 



PiT)^PqiTq), e{T)=eq{Tq) 

where H is to be regarded as the bulk pressure of d 
hydrodynamics which is defined by 

1 



H 



(3) 



Here 



({^ — is the four velocity difference be- 
tween the non-extensive parameter q > 1 and q — 1 
flows, 7(2:) = 5u'^{x)u^{x) = —^5uqfj_Su^ and Wq = 
Eq -\- Pq is the g-enthalpy. If T{x) and "f{x) satisfy 
Eq. ^ relations in the whole space-time region, one can 
always transform equation of perfect g-hydrodynamics 
Eq. (H]) into the following equation of d-hydrodynamics 



e{T)i 



'[P(T)-f H] A^'^+2M^K'')+7r''H =0, (4) 



in which one recognizes the energy flow vector and 
the (symmetric and traceless) shear pressure tensor i:^^ , 



(5) 



2 



(where A^^ e 
and a<^6''> 



. . , ^ _ 2. 
= [i(A^A^ + A;^A^) - iA^''AA.]a'^6^ 
The d-hydrodynamics represented by Eq. ^ can be 
therefore regarded as a viscous counterpart of the per- 
fect g- hydrodynamics represented by Eq. ([T]). We call 
this relation the nonextensive/dissipative correspondence, 
(NexDC). With the bulk pressure Eq. ^ and NexDC re- 



lations Eq. ([2|) one gets the q-enthalpy 



y/[l+7]' 



with w = Ts ^ e + P being the usual enthalpy in BG 
statistics. Accordingly, the bulk pressure Eq. ([3]) can be 
written as 



n = wT, where 



1 7(7 + 2) 
" 3(7 + iP' 



(6) 



Finally, the NexDC leads to the following relations be- 
tween components of the dissipative tensor: 



:6n2. (7) 



We shall now address the central question of this 
Letter: is d-hydrodynamics ^ obtained in this way 
causal, as is the one discussed in 

0, i, S, E, El? To 

this end, let us first notice that, because we expect 
that g — 1 <C 1 our stationary state described 

by Eqs. ([T]) and ([4]) can be regarded as a near equi- 
librium state (with energy momentum tensor 7^^" = 
{eq+Pq)u^u'^—Pqg^''). Accordingly, the state with q ~ I, 
i.e., the one with no residual correlations between fluid 
elements and no intrinsic fluctuations present, is a true 
equilibrium state (with energy momentum tensor TJ^ = 
£{T)u^'u'' - P(r)A'''' and with the equilibrium distri- 
bution given by the usual Boltzmann distribution func- 
tion, fcq{x,p) — e:xjp[—p'^u^{x)/kBT{x)]). Both states 
are near to each other and Tf^ — TJi'^ + 6T'^'^, where 
gq-fi,^ ^ -llA^"' + 2W'^^'u'''> +Tr>"' . As an immediate con- 
sequence of this, the conservation of g-entropy assumed 
in the ideal g-hydrodynamics, [squj^];^ = 0, results in the 
production of the usual entropy in d- hydrodynamics: 

, du 



T[su%, 



dr 



(^ = M^AT.p is the proper time derivative acting on a 
quantity X). Acting by d/dr on both sides of Eq. © 
and multiplying by 11, one obtains 



n 



dn dr ^ dw 
dr dr dr 



0. 



(9a) 



Similarly, acting d/dr on both sides of the first and the 
last terms in Eq. fT]) one gets two other identities: 



4"' 



dWf' 
dr 

~d^ 



n^ + nr^ 

dr dr 



0, 



^dn 

dr 



0. 



Finally, acting by (where V^X^ = AJ^X'^; 
sides of the middle term in Eq. ([7]), one gets: 



(9b) 

(9c) 
on both 

eO. (10) 



Following Q let us now introduce the usual thermody- 
namic (positive) coefficients, /3o,i,2i for the, respectively, 
scalar, vector and tensor dissipative contributions to the 
entropy current ([8|), and the viscous/heat coupling co- 
efficient a (all with dimension [GeV]"'' to ensure that 
dimensions of identities © and pH]) are the same as the 
dimension of r[sii''];^ in Eq. jH])). Multiplying ((9c)) . 
([9a|) and (fTO|) by the, respectively, /3o, Pi, (32 and Fa, and 
combining them together one gets the following identity: 

2 dW^" 

-(/3o -Pi- /32)n^ + 2arnyf,w^ 

1 dTr''" 



dr 

-(/3o+/3i)nF5^-t-/3ou;n^^=0 
dr dr 



(11) 



The g-entropy current is connected with su'^ and /T, 
1 T / „ W^'' 



T 



(12) 



Applying the four divergence to both sides of Eq. ([l3 
one gets (with cj)^ = 9^7/(1 + 7) + 9f_,Tq/Tq and $ = 



dw 
'd7 



(13) 



from which and Eq. ([5]) one has 



F^ = r(^w^^ - v.VK^) -n.^ + FVK''^. (w) 

dr dr 

(with M^.^ = V^W + and u^W^' = 0). Adding 

now Eq. pTjl to ([H]) and using Eq. one obtains 

Tisu^;^ = -nT-T4^^$^ + 7r'^'^*^, (15) 



(8) with T, $^ and defined as (/3^ = Po - Pi, PI = 



/3i+/3o, «* = ^+Pl): 

^-iP0-P2)^+U^, 

+ kP14>U - (2a - a*)FV^l^'^, 



__{Pl~_p^ dW^ 

^~ 3 dr 



i + a*nr)— ^ 

dr 



(16a) 

Pinr^,, 



(^L±M^w^^_«rv,.A_2^rv,n, (i6b) 



6 

_ P2 d-K^y 

PI 



6 



(1 - K)07r^^ +arV(^Wj,). 



(16c) 



The arbitrary constant k G (0, 1) appearing in Eqs. (I16a|) 
and (jl6cp is due to an ambiguity originating from the 
last NexDC relation in Eq. ([7]) which allows us to write 



3 



= kH^ + (1 - Hi)^T^"'■n^_,^/&. To ensure that [sm''];^ > 0, 
we assume now (following the standard 2nd order the- 
ory 0) a linear relationship between the thermodynamic 
fluxes n, W^^, irf^" and the corresponding thermodynamic 
forces, TjCJ^j^'piy, with the usual transport coefficients 
C, A, r], respectively. We have then 



2ri 



(17) 



and our g-hydrodynamics leads to the following consti- 
tutive equations of the corresponding d-hydrodynamics: 



n + rn— = 



dT 



dT 

■/?i*)c;r 



AT 



ni5l^<tA\^l^w^^W^, (18a) 

0?M„ 



(1 

Wtt VcrTT 



"*nr)^-/3*nr0^ 



j3* 

- K)Tj(l>n^„ 



/wnV^n,(18b) 



(18c) 



with the corresponding relaxation times rn = C(/3o ^ 
(32), Tw = |(/3r - Po)TXT and = |/32?7 and with 
the corresponding heat-viscous coupling lengths luw = 
C(2q; — a*)r, IwVl — 2Q;rAT, ^^/^r■K = aFAT and Z^vf = 
277ar. To ensure the positivity of relaxation times we 
must have /So > + /92 ■ Actually, one can also express 
coefficients /?o,i,2 and a in terms of the corresponding 
relaxation times and transport coefficients: 



_ rn Sttt 3 tw 
Q 7] 2 1 Aj 



P2 — , 



a ■ 



Iwtt 
TXT 



IttW 

2f^' 



3 Tw 



1 = 77 



2 TAT' 



(19) 



To summarize this part, under the conditions mentioned 
above, and with the help of the dissipative tensor re- 
lations, Eq. ([7]) and g-entropy conservation, one can 
obtain the corresponding constitutive equations of d- 
hydrodynamics, Eq. (|18p . They include the relaxation 
times and heat-viscous coupling lengths in a quite natural 
way. Since the original perfect q-hydrodynamics Eq. ([1]) 
does not contain any space-time scale, it is than natural 
that the NexDC conjecture does not introduce any defi- 
nite relaxation time or viscous-heat coupling length scale. 
However, as seen in Eq. (jl8[) . the d-hydrodynamics ob- 
tained from the g-hydrodynamics by using the NexDC re- 
lations, possesses the causal property of the correspond- 
ing d- hydrodynamics. Let us close noticing that this d- 
hydrodynamics predicts that 



21 



1 



XT ^ 
rrn) 



IttW , 



srr^) 



T}/S 



XT /a' 



(20) 
(21) 



As an example, let us now apply the NexDC to a q- 
hydrodynamics with the Bjorken type scaling initial con- 
ditions [131. In this case the equation of motion Eq. ([1]) 
for an q-ideal fluid has a very simple form: 



d£,(r) _ £,(r)+P,(r) 
dr T 



(22) 



with four velocity field uj^ = (1,0,0,0) (we use metric 
9tiu = {gTT,gxx,gvy,9r,",r ) = (1,0,0, -r^), where r = 
(^2 - z2)i/2 and 77* = iln(t + z)/{t - z)). The g-energy 
density and q-pressure are defined, respectively, as Eq = 



Uq^/T^^Uqi, and Pq = 
momentum tensor \\l 



-lT'"'Aqf_,^, by using the g-energy 



the fq(x,p) given by Tsallis distribution. 



fqix,p)^ i-(i-g) 



kBTq{x) 



1/(1-9) 



(23) 



The corresponding EoS, Pq — Pq{eq), is obtained by ehm- 
inating the common parameter Tq from both the Pq(Tq) 
and eq{Tq). Fig. [T] shows the proper time evolution of 
Tq, T (notice that, because of Eq. ([2|), they are not in- 
dependent) and the bulk pressure 11 for the relativistic 
TT gas (with mass m =0.14 GeV) distributed according 
to Eq. ((23)1 . The g-hydrodynamical evolution starts at 
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r[fiii] 
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FIG. 



1.00- 



1: Proper time evolution of T^, T and 11 for 
1.20. The g-hydrodynamical evolution starts at tq = 1.0 fm. 
The initial temperatures, Tq{To) are chosen in such way as to 
have To = 0.2 GeV for all values of q. The doted line in (b) 
is obtained for a fluid composed of massless particles (in this 
case there is no g-dependence) . 



proper time tq = 1.0 fm with the g-tcmpcrature Tqiro) 
chosen such that Tq = T(ro) = 0.2 GeV. As seen in 
Fig. m although Tq depends on g, the T{t) does not 
(it scales with g). This is interesting because in the 
usual d- hydrodynamics cooling behavior (i.e., T) is af- 
fected by the viscous effects. However, in our case the 
g-dependence of Eq. ([22]) is only implicit and the corre- 
sponding EoS depends only very we akly on g in the high 
temperature region (cf.. Fig. 1 in 13[). On the other 



4 



hand, it depends on the mass of particles which composed 
our fluid. The bulk pressure 11 decreases monotonously 
with r and increases with the non-extensive parameter 

q- 

We close this example by discussing the corresponding 
Reynolds number [2l[ . Separating in Eq. ([^ dissi- 
pative and non-dissipative terms by using the Reynolds 
number Q we get 



de e + P (e + P) 

:r + = ^ 

dT T 



In g-hydrodynamics with the NexDC conjecture 



3t 



P 



dn 
d7 



(24) 



(25) 



Fig. [2] shows examples of proper time evolution of the 
inverse of Reynolds number obtained in Eq. ((25| . For 
comparison, the maximal values of R~^ quoted in [9] for 
To = 1 fm and Tq = 0.3 GeV are R'^ « 0.065 (for 7]/s = 
0.08 and at r « 2 fm) and R-'^ w 0.21 (for 77/s = 0.30 
and T « 4 fm). Our R~^ remains smaller than these 
values, nevertheless we observe a similar maximum of 
R~^ at finite proper time, r ^ 10 fm. Notice that one 
can rewrite Eq. ([M]) as de/dr = {R~^ — 1)w/t, which, in 
the case of R^^ > L may indicate the instability of our 
solution for large r [7|]. However, as one can see in Fig. [21 
our q-hydrodynamical model with q-scaling solution gives 
very small values of Rj^, Rj^ ^ 1- This, according to 
[13, [Hi , niay guarantee the stability of the corresponding 
d-hydrodynamics obtained. 



(b) r = 0.3 GeV 



0^ 





^=1.05 



10 

T[fm] 



FIG. 2: Proper time evolution of R^^ obtained from q- 
hydrodynamics for two different initial temperatures (panels 
(a) and (b)) and for different values of the parameter q. Solid 
lines are for, respectively, q = 1.05, 1.10 and 1.20, dashed 
lines are for q = 1.30 — 1.50. The same EoS has been used as 
in Fig. m 



To summarize: we have demonstrated that 
d-hydrodynamics obtained from the perfect q- 
hydrodynamics (by means of the NexDC conjecture) 



preserves its original causality. We believe then that 
q-hydrodynamics can serve as a phenomenological model 
of d-hydrodynamic with parameter q describing sum- 
marily non ideality of the hadronic fluid. We have also 
shown that the corresponding inverse of the Reynolds 
number, 1/Re, may be small enough to ensure a stable 
hydrodynamical evolution of the corresponding strongly 
interacting quark/hadronic matter. 
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